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Abstract 


The  non-uniform  (non-linear)  equilibria  of  the  classical  (short  circuit)  Pierce  diode  and 
the  extended  (series  RLC  external  circuit)  Pierce  diode  are  described  theoretically,  and 
explored  via  computer  simulation.  It  is  found  that  most  equilibria  are  correctly  predicted 
by  theory,  but  that  the  continuous  set  of  equilibria  of  the  classical  Pierce  diode  at  L  =  2  ■/ 
are  not  observed.  The  stability  characteristics  of  the  non-uniform  equilibria  are  also  worked 
out,  and  are  consistent  with  the  simulations. 

Introduction 


J 


The  classical  Pierce  diode  is  a  theoretical  model  which  was  introduced  by  J.  R.  Pierce  [1]  to 
predict  the  maximum  electron  current  which  could  be  passed  through  a  plasma  device  without 
instability.  It  consists  of  two  parallel  planes  (plates  or  grids)  and  a  cold  electron  beam  traveling 
between  them.  The  electrons  are  neutralized  by  a  background  population  of  infinitely  massive  ions. 
(In  Pierce’s  model,  this  background  was  stationary,  but  for  the  purposes  of  this  article  it  is  better  to 
think  of  them  as  co-moving.)  The  planes  are  connected  by  a  wire,  and  so  are  at  the  same  potential. 
The  linearised  behavior  of  this  model  has  been  studied  in  detail  [2,3]. 

The  extended  Pierce  diode  is  a  similar  device  which  has  a  passive  circuit  element  —  either  a 
capacitor,  resistor,  or  inductor  —  in  place  of  the  short  circuit  between  the  electrodes  (see  Fig.  1). 
This  device  is  interesting  as  a  second  approximation  to  real  bounded  plasma  systems.  The  linear 
behavior  of  this  device  has  been  worked  out  by  Kuhn  and  Horhager  [4],  and  verified  by  simulation 

M. 

Godfrey  [6]  has  analysed  the  stationary  equilibria  for  the  classical  Pierce  diode,  and  I  shall 
describe  these  first.  Next  the  equilibria  of  the  extended  Pierce  diode  will  be  investigated  by  expanding 
on  the  formulas  derived  by  Godfrey.  Finally,  particle  simulations  will  be  presented  which  test  and 
support  the  predictions  of  the  theory. 
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In  his  investigations  of  the  classical  Pierce  diode,  Godfrey  has  discovered  an  elegant  integral 
equation  formulation  of  the  Pierce  diode  with  the  single  requirement  that  the  velocity  of  the  beam 
be  a  single-valued  function  of  position.  This  integral  equation  formulation  is  easily  modified  to  cover 
the  extended  Pierce  diode,  and  is  also  ideally  suited  for  linearization  about  any  equilibrium. 


Equilibrium  Theory 

Although  Godfrey  gives  the  equilibrium  solution,  he  does  not  derive  it,  so  I  will  give  a  brief 
derivation  here.  The  equilibrium  equations  are 


p(x)v(z)  =  pov0 

dE  _  p—  po 

dx  ~  e0 


(1) 

(2) 

(3) 


with  boundary  conditions 


P(0)  =  Po 

®(0)  =  v0 

f  E(z)dz  =  -V{E{ 0)) 
Jo 


(4) 

(5) 

(6) 


where  p,  v,  E,  V,  q/m,  and  <o  represent  the  electron  charge  density,  electron  velocity,  electric 
field,  potential  drop  across  the  diode,  electron  charge- to- mass  ratio,  and  the  dielectric  constant  of 
the  vacuum.  The  first  three  of  these  are  functions  of  the  position  z.  The  potential  drop  V  may 
be  a  function  of  the  electric  field  because  of  the  external  circuit  element.  For  the  classical  Pierce 
diode,  V  =  0.  The  potential  drop  cannot  depend  on  the  time  derivatives  of  the  electric  field,  since 
the  problem  is  a  static  one,  but  this  implies  that  there  can  be  no  current  flowing  in  the  external 
circuit  (unless  one  considers  the  constant  injection  flux  of  electrons  to  be  balanced  by  an  externa) 
current  instead  of  an  equal  flux  of  massive  ions).  This  in  turn  implies  that  the  external  resistance  or 
inductance  will  have  no  effect  on  the  equilibria.  Only  the  capacitive  case  will  produce  new  equilibria. 


vCvWv-: 
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The  equations  can  be  made  dimensionless  by  renormalizing  as  follows. 


a2  _  gpol*  _  f UpolV 
~  e0mv$  \  v0  / 


o’  =  £ 


*'  =  £ 


E'=  -$Le 


v  =  -4-v 

mvi 


Note  that  the  charge  of  the  electron  has  been  absorbed  into  the  electric  field,  so  that  the  electric 
field  is  of  the  opposite  sign  as  the  physical  electric  field. 

The  resulting  equations  (dropping  the  primes)  are 

fw  =  l  (7) 

dv 

v-  =  E  (8) 


with  boundary  conditions 


*0)  =  1 


v(0)  =  1 


f  E[z)dx  =  -V{E( 0)) 
Jo 


To  solve  these  equations  define  t  such  that 


and  1  s  0  at  i  =  0.  This  t  represents  the  time  it  took  a  fluid  element  (or  particle)  to  arrive  at 
position  z  from  its  time  of  injection.  Equations  8  and  9  become 

dv 

-  =  e  m 

^  s  os(l  -  v)  (14) 
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(Equation  7  was  used  to  derive  (14),  and  is  now  superfluous.)  Either  E  or  v  may  be  eliminated  to 
yield  a  harmonic  oscillator  equation.  The  solution,  taking  the  boundary  condition  (11)  into  account 


E  =  Eq  cos  at 


,  E0  .  , 

v  =  1  +  —  sin  or 
a 


x  ss  t+  — j(l  -  cos  at) 


The  constant  of  integration  Eq  represents  the  electric  field  at  the  injection  plane,  and  may  be  either 


positive  or  negative. 


A  new  condition  is  implicit  in  this  solution;  z  must  be  a  monotonically  increasing  function  of 


t,  so  that  da/dt  —  v  must  be  greater  than  zero.  Therefore  we  must  require  that 


1  +  — sin  at  >  0 
a 


for  0  <  t  <  T.  For  most  cases,  this  will  imply  the  condition  |£o|  <  <*•  Note  that  since  v  is  a  single¬ 


valued  function  of  t,  the  condition  (18)  is  enough  to  ensure  that  v  is  a  single- valued  function  of  x, 


which  is  necessary  for  the  fluid  approximation  to  be  valid. 


The  condition  that  v  be  single-valued  is  not  at  all  restrictive,  since  all  equilibria  of  interest  must 


itisfy  this  condition  anyway.  When  the  system  is  in  equilibrium,  the  total  energy  of  an  electron  is 


constant,  and  equal  to  l/2mv3  +  q<j>.  Thus  only  two  values  of  v  are  possible  at  any  given  position, 


implying  that  the  electron  trajectory  can  either  go  from  the  injection  plane  to  the  opposite  plane 


without  turning  around  at  all,  or  it  can  turn  around  once  and  return  to  the  injection  plane.  The 


case  in  which  the  electrons  are  turned  around  implies  a  large  potential  barrier,  which  cannot  be 


created  by  the  background  ions,  and  cannot  be  sustained  by  a  passive  external  circuit  element  in 


the  bee  of  the  large  ion  current.  Therefore,  only  the  case  in  which  v  is  a  single- valued  function  of  x 


is  of  interest. 


Equatioo  17  applied  at  x  =  1  gives  an  important  relation  between  the  time  T  that  a  particle  or 


fluid  element  takes  to  transit  the  system,  and  Eq  , 


T=  1-^(1 -cosoT) 


$ 

I 


« 


& 


> 

■v  .V 
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The  second  boundary  condition  (12)  can  now  be  applied.  Let  T  be  the  time  it  takes  a  fluid 
element  (or  particle)  to  transit  the  system.  The  new  boundary  condition  involves  the  potential  drop 


across  the  system, 


-V(E)  =  I*  E  dx 

Jo 

fT  dx 
=  /  E—dt 

Jo  dt 

=  J  Eoco&at  •  ^1  +  ~sinot^  dt 

-V(E)  =  ^  sin  aT  •  (l  +  ^  sin  aT^J 


Equilibria  without  a  Capacitor 

Let  us  first  consider  the  classical  Pierce  diode.  In  this  case,  V(E)  =  0,  so 


£o  =  0 


sin  aT  ~  0 


(i  +  |.i».r)  =  o 


The  first  case  («)  is  the  uniform  equilibrium,  and  the  third  case  (iii)  implies  that  v(l)  =  -1  which 
violates  the  condition  that  the  velocity  be  positive,  so  the  second  case  (»i)  is  the  interesting  one. 


The  solution  is  simply 


_ 

T  =  — 
a 


Putting  this  into  (19)  yields  two  results  depending  on  whether  n  is  even  or  odd.  If  n  is  even, 


and  Eg  i*  unconstrained  (except  for  condition  (18)).  If  n  is  odd, 

E0  =  |(o  -  nv)  (23) 

The  condition  that  the  velocity  be  positive  requires  (except  for  n  =  1  when  Eo  >  0)  that 
|£g|/a  <  1.  For  any  given  n  this  condition  limits  the  range  of  a  over  which  a  valid  solution 
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exists.  The  case  n  =  1  is  different  for  Eq  >  0,  since  sin  at  is  always  positive  in  this  case.  Thus 
for  n  ss  1,  there  is  no  upper  bound  on  Eq.  This  equilibrium  structure  can  be  diagrammed  as  in 
Godfrey  (see  Fig.  2).  The  stability  of  these  modes  will  not  be  dealt  with  systematically  yet,  but 
the  equilibria  with  Eo  <  0  are,  in  fact,  unstable  while  those  with  E0  >  0  are  stable.  If  the  initial 
value  of  the  electric  field  at  the  injection  plane  is  less  slightly  less  negative  than  the  value  Eq  for  an 
unstable  equilibrium,  then  the  system  will  decay  either  to  the  uniform  equilibrium  (if  the  uniform 
equilibrium  is  stable),  or  to  an  oscillating  state.  If  Eo  is  more  negative  than  the  equilibrium  value, 
then  the  growth  of  the  unstable  mode  requires  that  a  virtual  cathode  form. 

Equilibria  with  a  Capacitor 


When  the  externa]  circuit  contains  a  capacitor,  the  equations  become  more  complicated.  First, 
let  us  derive  the  proper  expression  for  ^(£(0))  (which  can  also  be  written  as  V(E0)).  The  voltage 


across  the  capacitor  is 


where  the  sign  of  Q  is  chosen  as  in  Fig.  1.  The  total  surface  charge  on  the  injection  plane  must  be 
<r  =  Q/A  (where  A  is  the  area  of  the  injection  plane)  if  the  unperturbed  state  Eo  —  V  =  0  is  to  be 
accessible.  This  surface  charge  gives  rise  to  the  electric  field  at  the  injection  plane,  Eo,  such  that 


Putting  all  this  together, 

V(Eo)  =  (^E0  (24) 

Converting  to  normalised  units  gives 


V'(EJ)=^ 


Where  Co  is  the  vacuum  capacitance  between  the  injection  and  collection  planes.  Thus,  defining 


yields  the  dimensionless  equation  (again  dropping  the  primes) 


(25) 


Equation  20  now  becomes 
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This  equation  can  be  solved  parametrically  by  setting  a  phase  variable  4>  =  aT,  and  combining  it 
with  (19).  The  result  is 

_  ^  •  ^«ina  4>  -  sin^  •  (1  -  cos#) 
a  ~  \  sin2  4>  +  ^(1  -  cos#) 


(28) 


Figure  3  shows  several  diagrams  like  Fig.  2  which  give  the  equilibrium  values  of  Eq  as  a  function 
of  a  for  various  values  of  C.  For  finite  values  of  the  capacitor,  it  can  be  seen  that  the  modes  which 
were  at  a  =  2nx  merge  with  the  modes  which  are  at  the  next  smallest  value  of  or.  For  large  values 
of  C  and  small  values  of  a,  this  merging  occurs  at  values  of  Eq  which  exceed  the  constraint  (18), 
and  so  are  not  seen  in  the  graphs.  When  the  graphs  are  extended  to  unphysical  values  of  Eq  (Fig. 
4),  it  is  seen  that  the  curve  doubles  back  in  order  to  connect  these  modes. 


Some  interesting  points  can  be  shown  analytically;  for  instance,  the  points  at  which  there  are 
bifurcations  with  the  uniform  equilibrium  (i.e.,  when  £o  — ►  0),  (26)  reduces  to 

sina  =  (27) 

(since  Eq  =  0  implies  T  =  1).  These  are  the  values  of  a  at  which  the  uniform  equilibrium  goes  from 
stable  to  unstable  behavior.  Note  that  this  equation  has  only  a  finite  number  of  roots,  implying  that 
beyond  a  certain  value  of  a  (for  a  given  C),  there  are  no  more  bifurcations. 


The  different  branches  of  the  curves  in  Fig.  3  correspond  to  ranges  of  #  which  are  bounded 
below  by  an  odd  multiple  of  r,  and  above  by  an  even  multiple  of  t.  Other  ranges  of  <f>  yield  a  very 
negative  £o.  As  might  be  expected,  as  C  — »  oo,  the  sloping  part  of  each  branch  corresponds  to  <j> 
near  to  the  odd  multiple  of  r,  and  the  part  of  the  branch  which  is  approaching  vertical  corresponds 
to  #  near  to  the  even  multiple  of  r.  It  is  useful  to  note  that  at  #  =  r,  the  limit  of  E/a  is  —s/2 
regardless  of  the  value  of  C.  While  this  solution  at  a  =  0  is  of  no  interest  physically,  since  it  violates 
condition  (18),  it  does  constrain  the  first  branch  of  the  curve. 


Simulation  Results 


The  particle  simulation  code  PD  W1  [7]  was  used  to  simulate  the  expected  equilibria.  As  with  the 
theory,  the  simulation  results  will  be  broken  up  into  classical  (short  circuit)  and  extended  (capacitive 
external  circuit).  The  simulation  parameters  are  shown  in  Table  I. 
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Classical  Pierce  Diode 

The  easiest  equilibria  to  simulate  are  the  stable  ones.  To  simulate  these,  the  uniform  equilibrium 
was  given  a  slight  perturbation,  and  the  equilibria  came  about  naturally.  Figures  5-7  show  the  phase 
space  plots  after  equilibrium  has  been  reached,  and  time  histories  of  the  electric  field  at  *  =  0  as 
the  simulation  approaches  equilibrium  for  a  equal  to  3x/2,  7t/2,  and  11t/2.  These  equilibria  were 
observed  by  Crystal  and  Kuhn  [3],  although  they  did  not  compare  them  with  theory.  While  in  all 
three  cases  the  electric  field  at  x  =  0  settles  down  rapidly  to  the  average  value  predicted  by  (23), 
the  results  are  not  without  surprises.  The  electric  field  at  x  =  0  for  both  a  =  7t/2  and  a  =  llx/2 
shows  oscillations  about  the  equilibrium  value  which  are  only  weakly  damped.  These  oscillations 
will  be  examined  shortly,  when  the  general  issue  of  stability  will  be  addressed. 

The  stable  equilibrium  with  n  =  1  extends  to  all  values  of  a  >  x  —  2.  This  can  also  be 
simulated  in  regions  where  other  equilibria  are  preferred,  such  as  a  =  5x/2.  To  do  this,  it  is  not 
enough  to  perturb  the  uniform  equilibrium,  since  this  perturbation  will  not  grow  to  the  desired 
equilibrium.  Instead,  the  electric  field  at  the  injection  plane  Eq  is  held  fixed  at  the  predicted 
equilibrium  value,  overriding  the  circuit  condition  for  one  transit  time.  This  constraint  should 
force  the  desired  equilibrium.  After  a  transit  time  (or  more),  the  constraint  can  be  removed,  and 
the  self-consistent  circuit  condition  reinstated.  The  result  for  a  =  5x/2  is  shown  in  Fig.  8.  There 
are  no  oscillations  about  the  equilibrium  value  of  Eq.  Table  II  summarizes  the  results  for  the  stable 
equilibria.  The  values  for  a  =  3x/2  and  a  —  5x/2  are  extraordinarily  accurate,  but  the  values  for 
a  =  7x/2  and  a  =  llx/2  seem  to  deviate  significantly.  Halving  the  grid  spacing  and  the  time  step 
reduced  the  error  for  a  =  llr/2  from  5%  to  1.5%,  and  Richardson  extrapolation  assuming  second 
order  accuracy  reduces  this  to  0.5%,  so  it  appears  that  the  deviation  is  of  numerical  rather  than 
physical  origin. 

The  unstable  equilibria  are  rather  difficult  to  simulate.  One  approach  is  to  again  fix  Eq,  over¬ 
riding  the  circuit  constraint,  at  a  value  either  above  or  below  the  predicted  equilibrium  value  of  Eq 
for  a  transit  time,  then  restore  the  circuit  constraint  and  see  whether  the  simulation  moves  away 
from  the  predicted  equilibrium.  This  was  done  for  the  a  =  r/2  equilibrium.  The  predicted  value 
for  Eo  in  equilibrium  is  Eq  =  -x3/8  »  -1.23,  so  two  simulations  were  run  at  initial  values  of  Eo 
at  —1.2  and  —1.3  (recall  that  the  sign  of  Eo  is  different  in  the  theory  from  the  simulations  with 
negatively  charged  electrons).  The  results  of  the  simulation  starting  at  Eq  =  -1.2  are  shown  in  Fig. 
9.  The  system  quickly  settles  to  the  stable  uniform  equilibrium  (Eq  =  0,  v  =  1).  When  Eq  =  -1.3, 
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however,  a  virtual  cathode  forms  (see  Fig.10).  If  the  duration  of  the  simulation  is  increased,  it  is 
seen  that  a  very  regular,  though  non-sinusoidal,  virtual  cathode  oscillation  has  set  in  (see  Fig.  11). 
Recall  that  this  oscillation  occurs  at  a  value  of  a  which  is  only  half  of  the  value  at  which  the  Pierce 
diode  first  becomes  linearly  unstable. 

The  unstable  equilibrium  at  a  =  5t/2  was  also  simulated.  The  results  are  shown  in  Fig.  12. 
Since  no  regular  virtual  cathode  oscillations  were  expected  (and  none  were  found),  Eq  was  again 
fixed  for  a  transit  time,  this  time  at  exactly  the  value  predicted  for  the  equilibrium.  There  are  two 
pieces  of  evidence  that  the  predicted  equilibrium  exists  at  very  nearly  the  predicted  value  of  E0. 
First,  the  potential  returns  almost  to  zero  at  x  —  1,  and  second,  when  released,  the  growth  away 
from  the  equilibrium  appears  to  be  exponential  rather  than  linear.  The  equilibrium  is  marred  by  the 
presence  of  two  trapped  electrons  (just  visible  near  x  =  .75  and  x  =  .85,  slightly  below  the  passing 
jlectrons).  According  to  continuum  theory,  these  trapped  particles  should  not  be  there,  but  the 
errors  inherent  in  simulation  allow  them  to  become  trapped  during  the  initial  transient  (the  diode  is 
initially  in  the  uniform  state).  If  these  particles  were  not  there,  the  agreement  might  be  even  better. 

The  equilibria  at  o  =  2 nr  are  especially  interesting  because  any  value  of  E0  which  does  not 
violate  (18)  should  be  a  valid  equilibrium.  Trouble  can  be  expected,  however,  since  the  dispersion 
relation  of  the  uniform  equilibrium  has  a  double  root  of  zero  at  a  =  2nx,  implying  a  secular 
instability  (linear  growth  with  time).  This  secular  instability  may  (and  the  section  on  stability 
will  show  shortly  that  it  does)  extend  to  the  non-uniform  equilibria.  This  secular  instability  is 
of  interest,  but  my  simulation  efforts  have  not  met  with  success.  Particle  simulation  at  a  =  2x 
produces  a  slowly  but  exponentially  growing  mode,  and  a  simulation  based  on  Godfrey’s  integral 
equation  formulation  produces  a  mode  which  oscillates  slowly  about  the  uniform  equilibrium  with 
large  amplitude.  The  reason  for  the  failure  of  these  simulations  seems  to  be  the  singular  nature  of 
the  a  =  2 nr  points.  For  values  of  a  near  these  points,  the  dominant  solution  of  the  dispersion 
relation  varies  as  9  ~  (2nx  -  a)1*2.  Thus,  a  very  small  numerical  error  can  produce  relatively 
large  deviations.  For  instance  an  error  which  alters  the  effective  value  of  a  by  1/1000,  may  create  a 
growth  rate  or  frequency  of  1/30,  which  is  quite  noticable  in  simulation. 

Equilibria  with  External  Capacitor 

To  test  the  predictions  for  the  equilibria  with  an  external  capacitor,  a  was  again  chosen  to  be 
half-integer  multiples  of  x  in  the  simulations.  The  results  for  the  stable  equilibria  are  in  Table  III. 


The  results  are  very  good  for  a  =  3x/2,  and  somewhat  less  good  for  a  =  7ir/2  As  was  demonstrated 
for  the  short  circuit  case,  this  less-accurate  agreement  is  most  likely  due  to  the  error  inherent  in  the 
simulation,  since  the  time  step  is  larger  relative  to  a  plasma  period  (i  e .,  is  larger) 

Interestingly,  the  equilibrium  for  a  =  7w/2  and  C  =  20  is  not  truly  stable  There  is  an  oscillatory 
mode  which  has  a  small  but  unmistakable  growth  rate  (see  Fig.  13)  This  growth  will  be  explained 
in  the  next  section. 

The  unstable  equilibrium  at  a  =  7x/2  and  C  =  10  was  also  simulated,  and  the  phase  space  plot 
is  shown  in  Fig.  14.  The  growth  rate  for  this  mode  is  quite  slow  relative  to  the  unstable  equilibria 
in  the  short  circuit  case. 

Stability  of  Non-uniform  Equilibria 

The  most  direct  method  of  analyzing  the  stability  of  an  equilibrium  is  to  compute  the  spectrum 
of  linear  perturbations  about  that  equilibrium.  Fortunately  this  is  possible  for  the  extended  Pierce 
diode.  One  simple  method  of  doing  this  is  to  extend  the  set  of  integral  equations  discovered  by 
Godfrey  for  the  classical  Pierce  diode.  Godfrey’s  integral  formulation  of  the  Pierce  diode,  extended 
to  the  case  with  an  external  circuit,  is  comprised  of  the  two  integro-difTerentia]  equations 

T(t)-l  =  -i  f  £(r)  sin  a(t  -  r)  dr  (28) 

a  Jt_T 

and 

(iS?r  +  b  +  0  E ^  =  T^1  “  ~a  Jt  TE^t  ~  r),,no(t  “  T^dT  <29) 

where  x  and  t  represent  position  and  time  respectively,  E  represents  the  e'  *  field  at  the  injection 
plane,  T  represents  the  transit  time  of  the  electron  just  leaving  the  system,  a  is  the  classical  Pierce 
parameter,  and  R,  L,  and  C  represent  the  external  resistance,  inductance,  and  capacitance  All  these 
quantities  have  been  normalized  (see  Appendix  A  for  a  complete  description  of  the  normalization 
factors  and  a  derivation  of  (28)  and  (29)). 

These  equations  are  idealy  suited  to  linearization  about  non-uniform  equilibria. 

Linearization  About  an  Equilibrium 

Assume  that  E  =  Eo,  T  =  To  is  a  solution  of  (28)  and  (29)  for  some  given  a  and  C,  then  (28) 
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and  (29)  can  be  linearized. 


If*  1  f*~To 

ST  = - /  6  E(r)  am  a(t  —  r)  dr - /  Eq  sin  a(t  —  r)  dr 

a  Jt-T  a  Jt-T 

«—  —  /  6E(r)  sma(t  —  t) dr  — —aiaaTo  ■  ST 

a  Jt-To  a 

^1  +  —  sin  aTo^j  ST  =  —  —  J  6  E(r)  sin  a(t  —  r)  dr 


The  equation  for  E(t)  can  be  similarly  linearized. 


{Lw+Rz+b+1)tE=  -°,T°ST  -  y(ST)i 

-  a  f  SE(r)(t  -  r)  sin  cr(t  -  r)  dr 

Jt-T 

rt- To 

—  ct  Eo(t  —  r)  sin  a(t  —  r)  dr 

Jt-T 

ss  —a7To  ^1  +  —  sin  q7o^  ST  —  SE(r)(t  —  r)  sin  a(t  —  r)  dr 


Substituting  in  (30), 

+  l^+l'jSE  ~a  J*  6E{r)(To  -  t  +  r)sina(t  -  r)dr 

=  a  /  *  6E(t  -  t')(T0  -  t1)  sin  at'  dt1  (31) 

Jo 

Both  sides  of  (31)  are  linear  in  6E,  so  a  solution  of  the  form  SE  —  exp  (dt)  can  be  sought.  When 
this  form  for  SE  is  substituted  into  (31),  the  result  is 

LO7  +  R0  +  ^  +  1  =  o  jfT#  e-#,'(T0  -  t') sint'dt' 

=  "  (ar|ga)»  12q*  C1  -  e"To#  «*qT°)  +  («’  -  ^)*inor0]  (32) 

Note  that  this  equation  does  not  contain  Eq. 

This  formula  (without  the  external  circuit  elements)  was  first  derived  by  Godfrey.  It  can  now 
be  applied  to  some  special  cases  of  interest. 

Uniform  Equilibrium 

In  the  uniform  (E  =  0  for  all  z)  equilibrium  To  =  1,  so 

LO7  +  R0+  —  +  aJ  +  0j  +  (qJ  +  02^  [2o5  (l  -  e~*  cos  a)  +(oJ  -  0J)sina]  =  0  (33) 
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This  is  the  dispersion  relation  found  by  Kuhn  and  Horhager  [4],  which  has  been  verified  by  simula¬ 
tions  [5]. 


Non-uniform  Equilibrium  with  Short  Circuit 


As  Godfrey  showed,  the  non-uniform  equilibria  when  the  external  circuit  is  a  short  (classical 


Pierce  diode)  are  of  two  classes  depending  on  whether  an  integer  parameter  n  is  odd  or  even.  When 
n  is  even,  the  equilibrium  is  given  by  o  =  nr,  To  =  1,  and  E0  with  any  value  as  long  as  -o  <  E0  <  a 


(as  per  Equation  22).  Putting  these  values  into  (32)  with  L  =  R  =  1/C  =  0  gives 


02  2  nJrJd  ,,  _#N  „ 

n’r’  +  9*  +  (n’x*  +  **)*  1  6 


Note  that  the  same  dispersion  equation  applies  for  all  values  of  £q.  The  dominant  root  of  this 


equation  is  9  —  0,  which  is  a  double  root,  implying  that  a  secular  instability  is  possible.  As  was 


mentioned,  this  secular  instability  was  not  observed  in  simulations,  and  the  singular  behavior  near 


a  =  nr  with  n  even  seems  to  be  responsible  for  the  poor  simulation  results. 


When  n  is  odd,  T0  =  nr/a,  and  £q  =  ^a(a  —  nr).  Substituting  these  values  into  (32)  gives 


Figures  15-17  show  the  dominant  solutions  of  this  dispersion  relation  as  a  function  of  a  for  the 
regions  of  interest  (nr  —  2  <  a  <  nr  +  2)  for  n  =  1, 3, 5.  For  nr  -  2  <  a  <  nr,  the  dominant  mode 


is  purely  growing  in  all  cases.  For  nr  <  a  <  nr  +  2,  the  dominant  mode  is  always  purely  damped 


for  n  =  1,  but  for  n  >  1  a  second  mode,  which  is  damped  and  oscillatory,  appears,  and  seems  to 


touch  Red  =  0.  This  is  not  an  illusion,  and  the  values  of  a  at  which  the  real  part  of  the  eigenvalue 


0  is  sero  can  be  found  by  setting  0  =  *w  and  setting  the  real  and  imaginary  parts  of  (35)  to  sero. 


The  result  is  a  Diophantine  equation, 


a  nm- 

- n  =  1 - 2 

r  n1  -m* 


in  which  m  must  be  odd.  Only  values  of  a/r  -  n  between  — 2/r  and  2/r  are  of  interest.  In  this 


region,  one  quickly  finds  that 


is  a  necessary  condition.  From  this  it  is  plain  that  for  n  s  1,  no  such  solutions  are  expected,  and 
for  n  =  3  and  n  =  5,  one  such  solution  is  expected.  Two  solution  will  not  appear  until  n  =  15. 


,  >  VV  V 


The  simulations  from  the  previous  article  agree  with  the  predictions  of  this  dispersion  relation 
as  well  as  can  be  determined  from  the  graphs  (the  eigenvalues  were  not  computed  numerically  from 
the  simulations). 

The  case  in  which  0  is  purely  imaginary  is  of  special  interest,  since  any  growth  or  decay  must 
be  determined  by  non-linear  effects.  Figure  18  shows  the  result  of  simulating  the  marginally  stable 
case  at  a  —  (3  4-  3/8)r.  The  simulation  moves  toward  equilibrium  and  quickly  reaches  a  steady 
oscillation  of  small  amplitude  about  the  equilibrium. 

Equilibria  with  an  External  Capacitor 

The  equilibria  with  external  capacitance  were  worked  out  earlier,  such  that  a,  To,  and  Eo  are 
all  functions  of  ^  =  a  To.  These  can  be  plugged  into  (32),  and  (32)  can  then  be  solved  numerically. 
Figures  19-21  show  the  eigenvalues,  9,  as  functions  of  a  for  several  different  values  of  C.  The  range 
of  a  in  the  plots  is  again  from  nr  —  2  to  nr  +  2,  although  this  is  no  longer  quite  the  right  range 
physically.  Note  that  the  oscillatory  modes  which  were  stable  in  the  short  circuit  case  have  become 
unstable  for  some  values  of  a  with  the  addition  of  an  external  capacitor,  counter  to  the  stabilizing 
influence  of  the  capacitor  on  the  modes  of  the  uniform  equilibrium.  These  theoretical  results  agree 
with  the  simulation  results  obtained  earlier. 

The  newly  unstable  modes  are  the  most  interesting  ones.  A  natural  question  is  whether  they 
saturate  at  some  finite  amplitude  due  to  non-linear  effects.  Figure  22  shows  that  for  C  =  20,  the 
a  =  7/2 x  case,  at  least,  does  not.  The  initial  transient  brings  the  system  close  to  the  (unstable) 
equilibrium  (at  t  ~  10),  but  then  the  instability  takes  over,  and  grows  until  finally  electrons  are 
turned  back  to  the  injection  plane.  Since  their  charge  remains  on  the  left  side  of  the  system  (and  the 
left  side  of  the  external  capacitor),  the  system  is  no  longer  in  a  state  which  allows  the  unperturbed 
(Eo  —  0,  V  =  0)  equilibrium,  and  so  the  resulting  equilibrium  (shown  at  t  =  64)  has  not  been 
included  in  my  analysis  of  equilibria.  (Such  equilibria  could  be  worked  out,  but  require  the  addition 
of  another  parameter,  specifically,  the  sum  of  the  charge  on  the  injection  plane  and  the  side  of  the 
capacitor  tied  to  it.) 

Once  a  certain  amount  of  charge  has  been  returned  to  the  emission  plane,  the  system  is  in  an 
equilibrium  which  has  a  stable  spectrum.  Interestingly,  the  potential  drop  across  the  system  in  this 
new  equilibrium  is  nearly  zero.  It  seems  unlikely  that  this  could  be  coincidence,  but  it  also  seems 


unlikely  that  a  grogs  process  such  as  virtual  cathode  formation  could  return  precisely  the  correct 
amount  of  charge  for  this  to  occur. 

Conclusion 

The  theory  for  non-uniform  (non-linear)  equilibria  of  the  Pierce  diode  with  an  external  capacitor 
instead  of  a  short  circuit  between  its  electrodes  was  worked  out,  and  simulations  were  performed 
which  verified  the  theory  for  both  the  capacitive  Pierce  diode,  and  the  classical  Pierce  diode. 

The  spectra  of  oscillations  about  non-uniform  equilibria  of  the  Pierce  diode  with  and  without 
an  external  capacitor  have  also  been  calculated  and  the  results  are  consistent  with  the  simulation 
results,  although  the  simulations  were  not  examined  to  high  precision. 

Of  particular  interest  is  the  observation  that  strong,  regular  virtual  cathode  cscillations  can 
occur  at  current  levels  which  are  much  smaller  than  the  current  level  at  which  the  Pierce  diode 
becomes  linearly  unstable  (simulations  were  performed  at  one  quarter  of  this  current  value).  As 
Godfrey  showed,  it  is  also  possible  to  exceed  the  current  level  at  which  the  Pierce  diode  becomes 
unstable  by  any  amount  desired,  if  the  proper  stable  non-uniform  equilibrium  can  be  attained. 

The  effect  of  the  external  capacitor  is  to  stabilize  the  linear  modes  somewhat,  and  also  to  increase 
the  current  value  at  which  virtual  cathode  oscillations  can  be  excited.  The  external  capacitance 
has  the  detrimental  effect,  however,  of  limiting  the  maximum  current  which  can  be  carried  by  the 
non-linear  modes.  The  most  useful  mode  (n  =  1)  ceases  to  exist  for  values  of  a  greater  than  some 
maximum,  and  the  other  modes  tend  to  become  unstable. 
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APPENDIX  A:  Derivation  of  Integral  Equations 


Renormalization  of  Equations 

The  physical  equations  of  evolution  are 


!+ £<„)=„ 

(al) 

{rt+'”k)',=  mE 

(°2) 

dE  p~  po 
dx  to 

(a3) 

and 

V-L££.+Rdc>  +  « 
v-Li?+Rir+c 

(a4) 

with  boundary  conditions  of,  respectively, 

fi(x  =  0)  =  po 

(a5) 

v(x  =  0)  =  Vo 

(a6) 

1) 

o 

II 

H 

S' 

(°7) 

['  Edx  =  -V 

Jo 

(a8) 

The  initial  conditions  must  also  be  specified.  The  independent  variables  are  x 

(position)  and  t 

(time).  The  dependent  variables  are  p  (charge  density),  v  (velocity),  E  (electric  field),  V  (voltage 

across  the  system),  and  Q  (charge  on  the  capacitor).  The  last  two  (V  and  Q)  depend  only  on  time. 

The  constants  in  the  equations  are  q/m  (electron  charge- to- mass  ratio),  to  (dielectric  constant  of 

the  vacuum),  po  (equilibrium  charge  density),  vo  (injection  velocity),  /  (length  of  the  system),  A 

(area  of  end  plates),  L  (external  inductance),  R  (external  resistance),  C  (external  capacitance),  Qq 

(initial  charge  on  capacitor),  7o  (initial  current  through  external  circuit).  The  condition  that  the 

electric  field  at  *  =  0  be  proportional  to  the  charge  on  the  capacitor  could  be  relaxed  to  allow  for  a 

constant  offset,  but  then  the  Q  =  E  —  0  state  would  not  be  an  allowed  state,  and  this  is  the  desired 
equilibrium. 

These  equations  can  be  renormalized  reducing  them  to  two  independent  variables,  three  inde¬ 
pendent  variables,  four  parameters  and  initial  conditions  as  follows.  First  define 
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Eo(t)  =  E{x  =  0,t) 


then  setting 


I  =  7 


>'  =  v4t 


/-£ 

Pa 


v'  =  ^- 
vo 


£'  =  -^y£ 
mi)n 


C=-j-C 

Aeo 


k  =  v-*£±r 


,  VqAcq  _ 
L  -“L 


yields  the  equations  (dropping  the  primes) 


(s+'£),,= 


37  =  »'('’ -l> 


£**m-(l£**i*h)* 


with  boundary  conditions 


=  0)  =  1 


v(x  =  0)  =  1 


Derivation  of  Integ^o- differential  Equations 


We  begin  by  taking  the  partial  derivative  with  respect  to  time  of  (all),  and  using  (a9)  to 
substitute  for  the  partial  derivative  of  p  with  respect  to  time.  The  result  is  the  familiar  (aside  from 
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the  appearance  of  a3)  equation 


d  (dE  ^  ,  \  , 

SVar  +  0,'’''J=’ 


This  equation  implies  that  the  quantity  inside  the  parentheses  is  a  function  only  of  time.  We  can 
therefore  take  its  value  for  any  value  of  z  to  be  its  value  at  z  =  0,  so  that 


*  +PV~  di  +a 


The  left-hand  side  of  this  equation  can  be  rewritten  as 


dE  dE  dE  (  dE\ 


=  it+av 


+  o3t> 


yielding 


j  dEo  2 
■  +  =  — r—  +  a 1 


This  equation  is  analogous  to  the  Llewellyn  equation  [8]  for  electron  beams  (v  =  Since 

v  =  dx/dt,  this  equation  can  be  integrated  once  with  respect  to  time.  Let  f0  be  the  time  at  which  a 
particle  (or  fluid  element)  at  time  1  and  position  x  was  emitted  at  z  =  0,  then 

£  t  -  Yi  +  «a*  =  ~  £o(to)  +  °l(t  ~  to) 

~  +  o3*  =  £b  +  oJ(‘-<o)  («15) 

The  boundary  conditions  on  this  equation  are 


z(t  =  t0)  =  0 


~  =1 

*  .... 


The  solution,  by  standard  methods  of  calculus,  is 


1  /‘  . 

=  t  -  to  +  —  /  Eo(r)  sin  a(t  -  r )  dr 
°  Jtm 
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Let  T  be  the  transit  time  of  the  particles  just  leaving  the  system,  i.e.,  z  =  1  when  t  —  tp  =  T.  Note 
that  T  is  a  function  of  time.  The  first  of  the  two  integro-differential  equations  to  be  derived  is 


1 

T  —  l  = - /  Ep(r)aina(t  —  r)dr 

°  Jt-T 


So  far,  to  has  been  viewed  as  a  parameter  in  a  solution  by  characteristics,  but  a  more  powerful, 
and  now  necessary  view  is  that  we  have  made  a  change  in  independent  variable.  To  consider  to  to 
be  an  independent  variable,  and  z  a  function  of  it,  it  is  necessary  only  to  show  that  the  differential 
operator  we  have  been  denoting  as  d/dt  is  actually  the  partial  derivative  with  respect  to  t  with  to 
held  fixed. 

—  —  -  L  ^fo  d_  dtp  d 

dt  .  +  V  dx  ~  dt  ..  +  dt  .  dtp  .  +  *  dx  .  dtp  , 


a 

dtp 

i  a.  w 

d_ 

'  fop 

+  t  -r—  • 

«  d*  t 

dtp 

dtp 

d 

i  §■ 

dtp 

(.  dx 

,  dt0  , 

dx 

We  are  now  free  to  make  use  of  to  as  an  independent  variable. 

Using  this  change  of  variable,  we  can  integrate  the  electric  field  to  find  the  potential  drop  across 
the  system,  in  order  to  convert  (a!2)  into  an  integro>dif!erential  equation. 


£  £(z,  t)  dz  =  J'~T  E(tp,  t)g~dtp 


Note  that  it  is  necessary  to  assume  at  this  point  that  z  is  a  monotonic  function  of  t0.  Since  dx/dtp 
is  negative  for  small  values  of  Ep,  it  follows  that  the  condition 

dx  1 

—  jg-  =  1  +  —  £o(to)sina(f  -  to)  >  0  (<*19) 

be  true  for  all  t  —  T  <  tp  <  t.  This  condition  will  always  be  satisfied  if  |£  +  0|  <  a. 

Since  z  is  known  explicitly  as  a  function  of  (  and  to,  and 

E-~ 

E~  dt* 

the  integral  can  be  reduced  to  a  single  explicit  integral.  The  algebra  is  lengthy  but  straight-forward, 
and  will  not  be  reproduced  here.  The  only  difficult  points  come  when  one  must  use  the  identities 


J [*  dx  I*  dy  f(y)  =  J'  dy  J'  dy  /(») 

=  j  (V  -  <*)/(»)  dy 


J'  dx ^  dy /(*)/( f(x)dx j 


The  final  result  is 


/l  E(x,t)dx  =  T*)  +  a  f  E0(r)(t  -  r)  sin  a(t  -  r)  di- 

Jo  2  y»-r 

Plugging  this  into  (al2)  gives  the  desired  integro-differential  equation 

(i5^,  +  A57+^+1)'E°=  T^l~T^~a  jt  T£o(T)(t-r)sin a(t-r)dr 


'  j'  v  ’/  V  V  V  A  A  A  ^ 
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EQUILIBRIA  FOR  EXTENDED  PIERCE  DIODE 


c 

Eo  Theory 

Eo  Simulation 

Q 

20 

-3.1399 

-3.1398 

3»/2 

10 

-2.5306 

-2.5306 

5 

-0.4863 

-0.4885 

20 

-5.3689 

-5.336 

7*/2 

Table  III.  Values  of  the  electric  field  at  the  cathode  for  different  values  of  a  and  C. 


ions 


electrons 


electrons 


ions 


electrons 


(c) 

Fi*  *.  Non-uniform  equilibria  for  (a)  C  *  30,  (b)  C  ■  10,  and  (e)  C  *  5.  Note  that  in  (a),  the  first  line 
dosa  not  fo  oo  indefinitely,  bat  ends  as  shown. 


(a)  (b) 


Fig.  5.  Phase  space  at  end  of  run  and  history  of  £«  for  short  circuit  case  with  a  —  3w/2.  (Phase  space  plot 
is  actually  comprised  of  many  particles,  and  not  a  single  curve,  as  it  appears.) 
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Fig.  13.  Phase  space  plot  and  history  of  Eo  for  alxnast-stabie  mode  at  a  s  7x/2  with  C  =  20 
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